Developing a Theoretical Framework for Optofluidic 
Device Designing for System IdentiHcation in Systems 
Biology: the EGFR Study Case 

Filippo MENOLASCINA^ , Vitoantonio BEVILACQUA\ Caterina CIMINELLi', 
Stefania TOMMASI^ and Angelo PARADISO^ 
^Polytechnic of Bari, Via E. Orabona 4, 70125, Italy 
'^Clinical Experimental Oncology Laboratory, National Cancer Institute, Via S. Hanhemann 

10, 70126, Bari, Italy 

f.menolascina@ieee. org 

Abstract: Identification of dynamics underlying biochemical pathways of interest in 
oncology is a primary goal in current systems biology. Understanding structures and 
interactions that govern the evolution of such systems is believed to be a cornerstone 
in this research. Systems theory and systems identification theory are primary 
resources for this task since they both provide a self consistent framework for 
modelling and manipulating models of dynamical systems that are best suited for the 
problem under investigation. We address herein the issue of obtaining an informative 
dataset to be used as starting point for identification of EGFR pathway dynamics. 
In order to match experimental identifiability criteria we propose a theoretical 
fi-amework for input stimulus design based on dynamical properties of the system 
imder investigation. A feasible optofluidic design has been designed on the basis of 
the spectral properties of the driving inputs that maximize information content after 
the theoretical studies. 
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1. Introduction 

Our understanding of molecular basis of complex diseases is being dramatically changed by 
systems investigation supported by the most advanced tools developed by the scientific 

community. In particular cancer has greatly benefited by systems level approaches since its 
development and progression is believed to be one of those system trajectories that arise 
Irom abnormal working states. The work by Romberg and colleagues [5] pointed out the 
relevance of a systems biology approaches to the study of dynamics bringing to cancer. 
Epidermal Growth Factor Receptor pathway is one of those biochemical pathways that is 
believed to play a central role in cancer development. Even if comprehensive models of 
EGFR pathway still exist [6], many questions still remain open in terms of parameter 
identifiability and driving input design in systems biology. In the next paragraphs we give a 
brief but formal description the problem under investigation and we describe a theoretical 
framework thought to provide researchers with effective guidelines to optimal experimental 
design in systems biology. 



2. Methods 



Dynamic models of biological systems can be conveniently represented by a general 
nonlinear system E„in the state space form depending on a generally unknown parameter 
vector p: 
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Where we have defined x as an n-dimensional state variable in which we store variables 
that generally represent some form of energy in our system; it is a m-dimensional vector 
that collects smooth input functions the system is supplied with; y is the r-dimensional 
observed output and p is a -dimensional vector of unknown parameters. This notation also 
assumes that /(•). 5'(0 and /i(-) are polynomial or rational functions of their parameters. 

A priori global identifiability of the system-experiment model, is a condition 
guaranteeing theoretical uniqueness of solution to the problem of recovering the model 
parameters from input-output data. Experimental identifiability, on the other hand, is 
strictly tied to the feasibility of a general identification process that starts from a pre- 
assigned set of input-output data. This last criterion can be seen as a primary key to develop 
guidelines for optimal experimental design since obtaining informative data. Given a quasi- 
stationary data set Z°° we say that Z°° is informative enough with respect to the model M* 
(the model to be identified) if, for any two models Wi{q) and W2(c[) in the set E\Wi{q) — 
W2{q)z{t)Y = which implies that W^^e'"^) = VKi2(e^'^) for almost all od. On the other 
hand a quasi-stationary data set Z°° is informative if it is informative enough with respect to 
the model set £*consisting of all linear, time invariant (LTI) models. The concept of 
informative dataset is tightly related to concepts of persistently exciting inputs. This can be 
seen easily by observing that a quasi-stationary dataset Z°° is informative if the spectrum 
matrix for z(t) = [u{f) y(t)]^ is strictly positive definite for almost all o). In fact, if we 
denote W-^^q) — Wiiq) = yV{q) applying a well known theorem on signal filtration we can 
write: /^J W{e^'^ j'&zC^J^) W^{e~>'^ ) doo = 0, with 4>x(<»^) spectrum of the signal x(f) and 
rOu(a)) Ouy(a))l 

where Oz((o) = ^ ^ ^ ^ ( \ ■ ^^^'■^^ ^z(t»i) is positive definite, this implies that 

V/(e^'^ ) = almost everjwhere that proves the previous statement. Moreover we can 
observe that, given the Oz(o)), for the Shur's Lemma we can assure algorithm convergence 
only if Ou(o)) > and Ou((o) - (o))Oyy (a))(&yu (oo) > 0. Evidently the only block of 
this array we have control on is the one representing the spectrum of input signal which 
directly depends on dynamical properties of the driving input signal we design. It is 
therefore convenient to introduce the concept of persistently exciting signal of order n for a 
quasi-stationary stimulus {u{f)}\ we say that a similar signal, with spectrum <I>u(co) is 
persistently exciting of order n is, for all filters of the form M^(q) = miq~^ + — h m^q'^ 
the relation IM^Ce^'^ )|^Ou((jl>) = implies that M^(e^'^) = 0. Evidently the function 
M^{z) Mji(^z~^) can have n — 1 different zeros on the unit circle (since one zero is always 
at the origin) taking symmetry into account. Hence u(t) is persistently exciting of order n if 
(.^) is different from zero on at least n points in the interval -n < o) < n. This is a 
direct consequence of the definition. Signals that show such properties have been 
investigated and include Pseudo Random Binary Signal, Generalized Binary Noise (or 
Random Binary Signal), Sum of Sines and Filtered Noise. Implementing such signals is 
quite simple from a computational point of view; however obtaining realizations of signals 



with such properties is an active area of research in current microfluidics. Developing 
geometries that satisfy physical conditions for the generation of signals compliant with the 
specifications imposed by the theoretical results we have shown is not a trivial task. Several 
alternative solutions have been proposed for signal modulation in microfluidic channels 
[1][2] [3]; they are based on diverse physical principles like boundary diffusion controlled 
by relative velocity (like in H filters), by exciting cells with diverse laminar fiows that 
affect different parts of the cell etc. Here we propose a method based on the theory of 
Kelvin-Helmholtz discontinuities (the idea is presented in Fig. 1). The effect results from 
velocity shears between two fiuid. Any time there is a non-zero curvature, the fiow of one 
fluid around another will lead to a slight centrifugal force which in turn leads to a change in 
pressure thereby amplifying the ripple. The most familiar example of this is wind blowing 
over calm water. Tiny dimples in the smooth surface will quickly be amplified to small 
waves and finally to frothing white-caps. Any sort of surface tension will hinder KH 
instabilities. If there is some restoring force Tg, the instability will arise if h.v^ > 
2(pi+pi)^y,^^^^ — On the other hand the stability requirement is met if bv^ < 2vj^ = 
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Moreover we can control the frequency and wave number of the instabilities by 



analyzing the dispersion equation that states that ^ = p+p^ 



p-p 



pp {U-U )2 



with 



k p+p' {p+p')^ 

p and p densities of the first and second fluid and U and U velocities of the flrst and 
second fluid respectively. 

3. Results and Discussion 

According to all these speciflcation we designed a microfluidic chip to be implemented 
using the polydimethilsiloxane (PDMS) technology. The mixing section of this chip is 
presented in Fig 2. 
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Fig. 1. Fluid behavior in the channel. Cells Fig. 2 Mixing section of the developed chip, 
are stimulated by a modulated signal 

As it can be seen in the final part of this design it thought to let fiuid converge in the 
common channel in which cells are anchored. Over this channel a here a fiber-coupled 
laser-diode module together with a photodetector allow to trace the presence and position of 
the cells; this ensures the . Computational fiuid dynamics studies on this geometry carried 
out using ANSYS CFD showed acceptable behaviors even at low Reynold regimes. Tests 
carried out using Potterswheel [4] MATLAB package to test theoretical results for the 
EGFR pathway (reference model by Kitano et al [6]) confirmed the expectations about 
information richness of the obtained I/O dataset and pushed the interest for further research 
in this field. 
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